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1. Introduction 

This Radcor symposium is devoted mainly to the high energy, LHC physics. However, low 
energy physics remains an active area of research. There are a few on-going low energy e + e~ 
experiments like the <I> factory DAONE (with the KLOE detector) the B factories PEP-II and KEK 
(with the BaBar detector at SLAC for which the accumulated data are still analysed and the Belle 
detector at KEK), the charm/r factory BEPC II (detector BES III). For them, to operate properly, a 
precise calculation of higher order corrections c t i wory for the process of Bhabha scattering (e + e~ — > 
e + e~) is necessary. Then the collider's luminosity / L tot dt = g m can be determined with high 
accuracy. Here N exp stands for a number of detected events, a is the theoretical cross section of 
the chosen reference process. If we know the luminosity accurately in a low energy region, then 
also the low energy hadron cross sections in e + e~ annihilation processes Chad = 7^ = ^theory j^- 
can be inspected {Nj w d is the number of measured hadronic events). What follows is that o t heory 
serves as the calibration parameter in measurements. 1 Two loop virtual QED 2 corrections to the 
Bhabha scattering are already known [^-|l^], however, so far their effects have not been applied to 
realistic physical situations (theoretical groups made comparisons only among virtual corrections, 
switching off the soft-photon cut-off parameter). For unresolved photons, fermion pairs, or hadrons 
at NNLO (real emission), 

e + e~ — > e + e~(y,YY), e + e~(e + e~), e + e~(f + f~), e + e~ (hadrons) , (1.1) 

we use the Fortran packages HELAC-PHEGAS (fermion pairs) |]17|— pCP, the MC event generator 



Ekhara [|21]-|24|] (pion pairs), and the Fortran program Bhagen-IPh-VAC (photons emission, 



M.Gunia, H.Czyz, unpublished, based on the generator Bhagen-IPh Q251). The complete results 



are compared with the BabaYaga [|26|-|28p MC generator in its recent version BabaYaga@NLO 



On the other hand, the muon pair production with real photon emission e + e~ — > jj, + jj,~Y is 
an important background and normalization reaction in the measurement of the pion form-factor: 
R exp = g(e+g-Z^+^-y) • We report on a progress in the determination of complete virtual corrections. 

2. NNLO Bhabha scattering at low energies 



In [30] it has been shown that for assumed realistic event selections the total NNLO massive 
corrections are relevant for precision luminosity measurements with 10~ 3 accuracy. Second, the 
maximum observed difference between BabaYaga@NLO and exact massive NNLO corrections 
is at the level of 0.07%. When cuts are varied the sum of the missing pieces can reach 0.1%, 
but for very tight acollinearity cuts only. The electron pair contribution is the largely dominant 
part of the correction. The muon pair and hadronic corrections are the next-to-important effects 
and quantitatively on the same grounds. The tau pair contribution is negligible for the energies of 



meson factories. Since the work [30] has already been reported at two conferences by M. Gunia 



("Matter To The Deepest", Ustroh 2011, [ETC]) and G. Montagna ("PHIPSI11", Novosibirsk 2011, 



in [g. 



'For more details on possible luminosity options at low energies, see the Introduction in [jip. 

2 Recently the dominant logarithmic 2-loop electroweak corrections, important at ILC energies, have been computed 
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[32]), we focus here on a part of results which has not been discussed there: main aspects of 
hadronic contributions and narrow resonances, for a complete discussion, see [|30|]. 

Specific problems are caused by narrow resonances like J/Y,\jf(2S), Typically, they are 

very narrow (at the keV level) and beam spread effetcs (at the MeV level) may cause that their final 
effects are diffused. This problem should be study. Here we show that their effects at fixed center 
of mass energy is dominating. 

Narrow resonances with mass M res and partial width T^J can be described approximately by 
the ansatz 



Rres(z) = ^M res rfj 8(Z- 



■Mi 



(2.1) 



Based on this, their contributions to the NNLO Bhabha process can be derived from the general for- 



mulae of QlSp. We discuss here as an example the contribution from the "rest" (genuine irreducible 



2-loop corrections, as defined in Q15|]). According to Eq. (87) of Jl5| ] it reads: 
9a 



dOres\ 

dQ. 



,2 T^e + e 
res 



K s M K 



+ ■ 



1 

res 



F 2 (M4)+F 3 (M4)ln 



1 



M 2 

res 



(2.2) 



Let us note that Eq. becomes invalid when the center of mass energy comes too close to 
the position of a resonance, i.e. if (s — M 2 es ) < F^/ M res . In the numerical examples, Table [l| this 
is not the case. 



Table 1: Soft+virtual NNLO contributions (J^l° es from narrow resonances (n.r.) defined by Eq. (2.2) for 



the Bhabha process with o/-Ebeam = 10~ 4 (in nb). The narrow resonance located closest to the center of 
mass energy of the given collider is included (first column, res) and excluded (second column, res'). The 
third column contains the Born cross section (Tg. 





U NNLO 
u rest,res 


—.NNLO 

rest, res' 


o B 


KLOE 1.020 


[all n.r.] 
-0.04538 


[n.r. without J/y/{lS)] 
-0.0096 


529.5 


BES 3.097 


[all n.r.] 
228.08 


[n.r. without J/yr(lS)] 
-0.0258 


14.75 


BES 3.650 


[all n.r.] 
-0.1907 


[n.r. without y/(2S)] 
-0.023668 


123.94 


BES 3.686 


[all n.r.] 
-62.537 


[n.r. without y/(2S)] 
-0.0254 


121.53 


BaBar 10.56 


[all n.r.] 
-0.0163 


[n.r. without T(4S)] 
-0.01438 


6.744 


Belle 10.58 


[all n.r.] 
0.04393 


[n.r. without T(4S)] 
-0.0137 


6.331 



In Table [T] we show numerical results based on Eq. (2.2). We can see that the contributions 
from narrow resonances dominate the NNLO Bhabha correction for BES running at J/y and y(2S) 
energies. For the remaining cases narrow resonances contribute below the per mille level when 
compared to the Born cross section Ob- For those cases results can be found in [30]. 
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We conclude that for experiments performed on top of a narrow resonance, this resonance can- 
not be treated as a mere correction and more detailed studies have to be performed. These should 
include the examination of finite width effects, beam spread effects, the estimation of NNNLO 
corrections, and the accuracy of the vacuum polarisation insertions in a close vicinity of these 
resonances. 

3. e + e — > n + /i y at low energies 

The e + e~ — > pL + pL~y process is an ideal benchmark process to test for massive tensor reduc- 
tions at the 1-loop level. It has (i) two different masses (and thus large difference of scales, up to 
7 orders in magnitude); (ii) quasi-collinear regions (due to small electron mass); (iii) small number 
of diagrams. 



+ 



b 




Figure 1: A set of diagrams needed for final state radiation (FSR) gauge invariance of interferences between 
tree diagrams (upper picture) and four and five point one-loop integrals (below). Here the diagrams are 
limited to FSR cases, the same property holds if initial state radiation (ISR) amplitudes are present. 

There are many tools which give a possibility to obtain automatically cross sections at one- 
loop level, however, there is still ongoing progress as new demands on precision of calculations 
require a more and more precise numerical analysis. Here we give an example of such a progress 
in the calculation of effects including 5-point functions. In Fig.[j] a class of diagrams is shown 
needed for building gauge invariant amplitudes. Let us define the relative accuracy 



max 



I^A^(<o P M tL) 
min9*(M< oop M t + ree ) 



(3.1) 



where i stands for (a), (b),(c) inFig.|I[ For such a test parameter, with the package LoopTools/FF 
[33, 34], 2.5 x 10 6 events have been generated. For Fortran double precision, an accuracy at the 
level A = 10~ 2 has been obtained, and for quadruple precision (with an Intel Fortran compiler), it 
is A = 10~ 12 . It is clear that using quadruple precision, a reasonable level of accuracy has been 
obtained. In Fig.|| a distribution of muon pairs in e + e~ — > /x+jU^y is shown at KLOE kinematics. 3 
Parameters as defined in Table || have been used. The results in FigsJ^ and |] estimate only the 

3 Similar plots for the same number of generated events but different distributions can be found in [B5I]). 
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Figure 2: Muon pair distributions including 5-point functions at KLOE. 2.5 x 10 6 events have been gener- 
ated. Looptools and FF packages have been used. 
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Figure 3: Muon pair distributions including 5-point functions at KLOE. 10 9 events have been generated. 
Looptools and FF packages have been used. 



virtual contributions. However, with an increasing number of generated events some instabilities 
appear, see Fig.||. Since version 2.2, LoopTools allows to choose in which decomposition five 
point integrals should be calculated. Implementations of both the Passarino-Veltman [ j36[ ] and the 
Denner/Dittmaier [|7|, 38] schemes exist. The instability problems are due to the appearance of 
inverse Gram determinants, and their solution is achieved with the alternative reduction scheme 
[39-41] and its implementation in the package P JFry. 4 



Ecms 






Q 2 




1.02 GeV 


0.02 GeV 


0°-15°, 165°-180° 


0.25-1.06 GeV 2 


50°-130° 



Table 2: Phase-space cuts for KLOE settings used in calculations. Q 2 is the invariant mass squared of the 
muon pair. 



Using these new reductions, in Fig.[| the same stability of distributions as in Figs.|2||3| is ob- 
tained. The results are completely stable and well controlled. E.g. the leading inverse Gram 
determinants \G^\ are eliminated in the reduction and small inverse Gram determinants \G^\ are 
avoided using asymptotic expansion. 



4 For the evaluation of the scalar integrals we use QCDloop/FF J4^]. For further references and some basic compar- 
431E34 
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Figure 4: Muon pair distributions including 5-point functions at KLOE. Bottom: absolute error estimate. 
Approximately 4 • 10 10 events has been generated. PJFry package has been used. Taken from Q]. 



4. Summary 

Low energy physics is an important field of activity. Let us mention the measurement of 
(g — 2)^ and its theoretical calculation, 5 but also low energy hadron physics, e.g. the determination 



of form factors. In this context, an important quantity is the pion form factor, see e.g. Q45|]. In or- 
der to describe it properly, experimental data are needed, and for that the process e + e~ — > /i+jU^y 
serves as normalization reaction. We have shown some details on progress in the calculation of 
this process at the NLO level, including 5-point diagrams. The progress has been possible due to 
a stable treatment of tensor reductions. Similarly, theoretical progress in calculations of NNLO 
effects in Bhabha scattering convenience us that the existing generator BabaYaga@NLO is a reli- 
able tool used for luminosity determination at meson factories. Here we have signaled only that 
some additional work on the understanding of narrow resonance contributions to precise Bhabha 
scattering studies would be welcome. 
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